New Vibrio cholerae sequences from Eastern and Southern Africa alter our understanding of regional cholera transmission.

Despite ongoing containment and vaccination efforts, cholera remains prevalent in many countries in sub-Saharan Africa. Part of the difficulty in containing cholera comes from our lack of understanding of how it circulates throughout the region. To better characterize regional transmission, we generated and analyzed 118 Vibrio cholerae genomes collected between 2007-2019 from five different countries in Southern and Eastern Africa. We showed that V. cholerae sequencing can be successful from a variety of sample types and filled in spatial and temporal gaps in our understanding of circulating lineages, including providing some of the first sequences from the 2018-2019 outbreaks in Uganda, Kenya, Tanzania, Zambia, and Malawi. Our results present a complex picture of cholera transmission in the region, with multiple lineages found to be co-circulating within several countries. We also find evidence that previously identified sporadic cases may be from larger, undersampled outbreaks, highlighting the need for careful examination of sampling biases and underscoring the need for continued and expanded cholera surveillance across the African continent.

Integrating dynamical modeling and phylogeographic inference to characterize global influenza circulation.

Global seasonal influenza circulation involves a complex interplay between local (seasonality, demography, host immunity) and global factors (international mobility) shaping recurrent epidemic patterns. No studies so far have reconciled the two spatial levels, evaluating the coupling between national epidemics, considering heterogeneous coverage of epidemiological and virological data, integrating different data sources. We propose a novel combined approach based on a dynamical model of global influenza spread (GLEAM), integrating high-resolution demographic and mobility data, and a generalized linear model of phylogeographic diffusion that accounts for time-varying migration rates. Seasonal migration fluxes across global macro-regions simulated with GLEAM are tested as phylogeographic predictors to provide model validation and calibration based on genetic data. Seasonal fluxes obtained with a specific transmissibility peak time and recurrent travel outperformed the raw air-transportation predictor, previously considered as optimal indicator of global influenza migration. Influenza A subtypes supported autumn-winter reproductive number as high as 2.25 and an average immunity duration of 2 years. Similar dynamics were preferred by influenza B lineages, with a lower autumn-winter reproductive number. Comparing simulated epidemic profiles against FluNet data offered comparatively limited resolution power. The multiscale approach enables model selection yielding a novel computational framework for describing global influenza dynamics at different scales - local transmission and national epidemics vs. international coupling through mobility and imported cases. Our findings have important implications to improve preparedness against seasonal influenza epidemics. The approach can be generalized to other epidemic contexts, such as emerging disease outbreaks to improve the flexibility and predictive power of modeling.

Synchrony of Bird Migration with Global Dispersal of Avian Influenza Reveals Exposed Bird Orders.

Highly pathogenic avian influenza virus (HPAIV) A H5, particularly clade 2.3.4.4, has caused worldwide outbreaks in domestic poultry, occasional spillover to humans, and increasing deaths of diverse species of wild birds since 2014. Wild bird migration is currently acknowledged as an important ecological process contributing to the global dispersal of HPAIV H5. However, this mechanism has not been quantified using bird movement data from different species, and the timing and location of exposure of different species is unclear. We sought to explore these questions through phylodynamic analyses based on empirical data of bird movement tracking and virus genome sequences of clade 2.3.4.4 and 2.3.2.1. First, we demonstrate that seasonal bird migration can explain salient features of the global dispersal of clade 2.3.4.4. Second, we detect synchrony between the seasonality of bird annual cycle phases and virus lineage movements. We reveal the differing exposed bird orders at geographical origins and destinations of HPAIV H5 clade 2.3.4.4 lineage movements, including relatively under-discussed orders. Our study provides a phylodynamic framework that links the bird movement ecology and genomic epidemiology of avian influenza; it highlights the importance of integrating bird behavior and life history in avian influenza studies.

Contemporary and historical human migration patterns shape hepatitis B virus diversity.

Infection by hepatitis B virus (HBV) is responsible for approximately 296 million chronic cases of hepatitis B, and roughly 880,000 deaths annually. The global burden of HBV is distributed unevenly, largely owing to the heterogeneous geographic distribution of its subtypes, each of which demonstrates different severity and responsiveness to antiviral therapy. It is therefore crucial to the global public health response to HBV that the spatiotemporal spread of each genotype is well characterized. In this study, we describe a collection of 133 newly sequenced HBV strains from recent African immigrants upon their arrival in Belgium. We incorporate these sequences-all of which we determine to come from genotypes A, D, and E-into a large-scale phylogeographic study with genomes sampled across the globe. We focus on investigating the spatio-temporal processes shaping the evolutionary history of the three genotypes we observe. We incorporate several recently published ancient HBV genomes for genotypes A and D to aid our analysis. We show that different spatio-temporal processes underlie the A, D, and E genotypes with the former two having originated in southeastern Asia, after which they spread across the world. The HBV E genotype is estimated to have originated in Africa, after which it spread to Europe and the Americas. Our results highlight the use of phylogeographic reconstruction as a tool to understand the recent spatiotemporal dynamics of HBV, and highlight the importance of supporting vulnerable populations in accordance with the needs presented by specific HBV genotypes.

Application of Phylodynamic Tools to Inform the Public Health Response to COVID-19: Qualitative Analysis of Expert Opinions.

BACKGROUND: In the wake of the SARS-CoV-2 pandemic, scientists have scrambled to collect and analyze SARS-CoV-2 genomic data to inform public health responses to COVID-19 in real time. Open source phylogenetic and data visualization platforms for monitoring SARS-CoV-2 genomic epidemiology have rapidly gained popularity for their ability to illuminate spatial-temporal transmission patterns worldwide. However, the utility of such tools to inform public health decision-making for COVID-19 in real time remains to be explored. OBJECTIVE: The aim of this study is to convene experts in public health, infectious diseases, virology, and bioinformatics-many of whom were actively engaged in the COVID-19 response-to discuss and report on the application of phylodynamic tools to inform pandemic responses. METHODS: In total, 4 focus groups (FGs) occurred between June 2020 and June 2021, covering both the pre- and postvariant strain emergence and vaccination eras of the ongoing COVID-19 crisis. Participants included national and international academic and government researchers, clinicians, public health practitioners, and other stakeholders recruited through purposive and convenience sampling by the study team. Open-ended questions were developed to prompt discussion. FGs I and II concentrated on phylodynamics for the public health practitioner, while FGs III and IV discussed the methodological nuances of phylodynamic inference. Two FGs per topic area to increase data saturation. An iterative, thematic qualitative framework was used for data analysis. RESULTS: We invited 41 experts to the FGs, and 23 (56%) agreed to participate. Across all the FG sessions, 15 (65%) of the participants were female, 17 (74%) were White, and 5 (22%) were Black. Participants were described as molecular epidemiologists (MEs; n=9, 39%), clinician-researchers (n=3, 13%), infectious disease experts (IDs; n=4, 17%), and public health professionals at the local (PHs; n=4, 17%), state (n=2, 9%), and federal (n=1, 4%) levels. They represented multiple countries in Europe, the United States, and the Caribbean. Nine major themes arose from the discussions: (1) translational/implementation science, (2) precision public health, (3) fundamental unknowns, (4) proper scientific communication, (5) methods of epidemiological investigation, (6) sampling bias, (7) interoperability standards, (8) academic/public health partnerships, and (9) resources. Collectively, participants felt that successful uptake of phylodynamic tools to inform the public health response relies on the strength of academic and public health partnerships. They called for interoperability standards in sequence data sharing, urged careful reporting to prevent misinterpretations, imagined that public health responses could be tailored to specific variants, and cited resource issues that would need to be addressed by policy makers in future outbreaks. CONCLUSIONS: This study is the first to detail the viewpoints of public health practitioners and molecular epidemiology experts on the use of viral genomic data to inform the response to the COVID-19 pandemic. The data gathered during this study provide important information from experts to help streamline the functionality and use of phylodynamic tools for pandemic responses.

Reconstruction of the origin and dispersal of the worldwide dominant Hepatitis B Virus subgenotype D1.

Hepatitis B is a potentially life-threatening liver infection caused by the hepatitis B virus (HBV). HBV-D1 is the dominant subgenotype in the Mediterranean basin, Eastern Europe, and Asia. However, little is currently known about its evolutionary history and spatio-temporal dynamics. We use Bayesian phylodynamic inference to investigate the temporal history of HBV-D1, for which we calibrate the molecular clock using ancient sequences, and reconstruct the viral global spatial dynamics based, for the first time, on full-length publicly available HBV-D1 genomes from a wide range of sampling dates. We pinpoint the origin of HBV subgenotype D1 before the current era (BCE) in Turkey/Anatolia. The spatial reconstructions reveal global viral transmission with a high degree of mixing. By combining modern-day and ancient sequences, we ensure sufficient temporal signal in HBV-D1 data to enable Bayesian phylodynamic inference using a molecular clock for time calibration. Our results shed light on the worldwide HBV-D1 epidemics and suggest that this originally Middle Eastern virus significantly affects more distant countries, such as those in mainland Europe.

SARS-CoV-2 introductions and early dynamics of the epidemic in Portugal.

Background: Genomic surveillance of SARS-CoV-2 in Portugal was rapidly implemented by the National Institute of Health in the early stages of the COVID-19 epidemic, in collaboration with more than 50 laboratories distributed nationwide. Methods: By applying recent phylodynamic models that allow integration of individual-based travel history, we reconstructed and characterized the spatio-temporal dynamics of SARS-CoV-2 introductions and early dissemination in Portugal. Results: We detected at least 277 independent SARS-CoV-2 introductions, mostly from European countries (namely the United Kingdom, Spain, France, Italy, and Switzerland), which were consistent with the countries with the highest connectivity with Portugal. Although most introductions were estimated to have occurred during early March 2020, it is likely that SARS-CoV-2 was silently circulating in Portugal throughout February, before the first cases were confirmed. Conclusions: Here we conclude that the earlier implementation of measures could have minimized the number of introductions and subsequent virus expansion in Portugal. This study lays the foundation for genomic epidemiology of SARS-CoV-2 in Portugal, and highlights the need for systematic and geographically-representative genomic surveillance.

Proposal for Human Respiratory Syncytial Virus Nomenclature below the Species Level.

Human respiratory syncytial virus (HRSV) is the leading viral cause of serious pediatric respiratory disease, and lifelong reinfections are common. Its 2 major subgroups, A and B, exhibit some antigenic variability, enabling HRSV to circulate annually. Globally, research has increased the number of HRSV genomic sequences available. To ensure accurate molecular epidemiology analyses, we propose a uniform nomenclature for HRSV-positive samples and isolates, and HRSV sequences, namely: HRSV/subgroup identifier/geographic identifier/unique sequence identifier/year of sampling. We also propose a template for submitting associated metadata. Universal nomenclature would help researchers retrieve and analyze sequence data to better understand the evolution of this virus.

Novel hepatitis B virus subgenotype A8 and quasi-subgenotype D12 in African-Belgian chronic carriers.

BACKGROUND: Hepatitis B virus (HBV) is a public health threatening virus and is classified into more than eight genotypes and more than forty subgenotypes. OBJECTIVES: To characterize and propose novel strains assigned as A8 and D12. METHODS: Four out of 133 HBV complete genome sequences, isolated from Belgian chronic carriers with African origin were phylogenetically analyzed. RESULTS: Phylogenetic analyses of HBV genotypes A and D strains exhibited separate clusters supported by significant bootstrap values. The two genotype A strains isolated from Congolese patients, and two genotype D strains isolated from Ghanaian carriers clustered separately from the other known subgenotypes A (A1-A6 and quasi-subgenotypes) and subgenotypes D (D1-D11). The mean inter-subgenotypic nucleotide divergence over the full-length genome sequence between the novel strains (A8 and D12) and A1-A7 and D1-D11 subgenotypes was higher than 4%. CONCLUSIONS: Phylogenetic analysis of the full-length HBV genome sequences revealed a novel subgenotype and quasi-subgenotype based on the nucleotide divergence and identification of novel amino acids motifs in different ORFs. We identified two strains of the novel subgenotype A8 and two strains of the novel quasi-subgenotype D12. Notably, the analysis demonstrated that the subgenotype A8 strains are a basal lineage that diverged before the other African subgenotypes A.

Evolution of highly pathogenic H7N3 avian influenza viruses in Mexico.

Highly pathogenic H7N3 influenza A viruses have persisted in poultry in Mexico since 2012, diversifying into multiple lineages that have spread to three Mexican states, as of 2016. The H7N3 viruses segregate into three distinct clades that are geographically structured. All 2016 viruses are resistant to adamantane antiviral drugs and have an extended 24-nucleotide insertion at the HA cleavage site that was acquired from host 28S ribosomal RNA.

Molecular epidemiology, phylogenetic analysis and genotype distribution of hepatitis B virus in Saudi Arabia: Predominance of genotype D1.

Despite the implementation of various vaccination programs, hepatitis B virus (HBV) poses a considerable health problem in Saudi Arabia. Insight on HBV evolutionary history in the region is limited. We performed a comprehensive epidemiological and phylogenetic reconstruction based on a large cohort of HBV infected patients. Three hundred and nineteen HBV-infected patients with different clinical manifestations, including inactive and active chronic carriers and patients with cirrhosis and hepatocellular carcinoma (HCC), were enrolled in this study. The full-length large S gene was amplified and sequenced. Phylogenetic analysis was performed to determine the genotype and subgenotypes of the isolates. Phylogenetic tree analysis revealed that genotype D is the most dominant genotype among patients. Moreover, this analysis identified two strains with genotype E isolated from active carriers. Detailed phylogenetic analyses confirmed the presence of four HBV D subgenotypes, D1 (93%, n = 296), D2 (0.02%, n = 5), D3 (0.003%, n = 1), and D4 (0.003%, n = 1). In addition, six genotype D strains were not assigned to any existing HBV D subgenotype. The large S gene of eight strains showed signatures of genotype recombination between the genotypes D and A and between D and E. Several strains harbored medically important point mutations at the protein level. Along with the dominance of the HBV genotype D, isolation of the E genotype and several recombinant strains from patients with Saudi Arabian origin is an essential result for decisions involving therapeutic measures for patients. Development of vaccines and detection of diagnostic escape mutations at antigenic epitopes on the HBsAg will be valuable to public health authorities. Furthermore, the diversity at the nucleotide and amino acid levels and different proportions of dN/dS at the PreS1, PreS2, and HBsAg reveal the selective pressure trend from inactive status towards advanced liver diseases.

15year fulminant hepatitis B follow-up in Belgium: Viral evolution and signature of demographic change.

Fulminant hepatitis among different clinical outcomes of hepatitis B virus infection is very rare and manifests high mortality rate, however it has not been investigated in Belgian inhabitants yet. In the frame of a retrospective study between 1995 and 2010, 80 serum samples (in some cases serial samples) archived in Biobank, were collected from 24 patients who had clinically developed fulminant infection of hepatitis B virus. In total, 33 hepatitis B virus (HBV) strains (31 full-length genome and 2 partial viral genes) of different HBV genotypes and subgenotypes including A2, B2, D1, D2, D3 and E, were amplified, sequenced and phylogenetically analyzed. HBV isolated strains from native and exotic patients were characterized by genome variations associated with viral invasiveness. Although several mutations at nucleotide and protein levels were detected, evolutionary analyses revealed a negative selective pressure over the viral genomes. This study revealed influence of immigration through a steady change in the viral epidemiological profile of the Belgian population.

Quantifying Next Generation Sequencing Sample Pre-Processing Bias in HIV-1 Complete Genome Sequencing.

Genetic analyses play a central role in infectious disease research. Massively parallelized "mechanical cloning" and sequencing technologies were quickly adopted by HIV researchers in order to broaden the understanding of the clinical importance of minor drug-resistant variants. These efforts have, however, remained largely limited to small genomic regions. The growing need to monitor multiple genome regions for drug resistance testing, as well as the obvious benefit for studying evolutionary and epidemic processes makes complete genome sequencing an important goal in viral research. In addition, a major drawback for NGS applications to RNA viruses is the need for large quantities of input DNA. Here, we use a generic overlapping amplicon-based near full-genome amplification protocol to compare low-input enzymatic fragmentation (Nextera™) with conventional mechanical shearing for Roche 454 sequencing. We find that the fragmentation method has only a modest impact on the characterization of the population composition and that for reliable results, the variation introduced at all steps of the procedure--from nucleic acid extraction to sequencing--should be taken into account, a finding that is also relevant for NGS technologies that are now more commonly used. Furthermore, by applying our protocol to deep sequence a number of pre-therapy plasma and PBMC samples, we illustrate the potential benefits of a near complete genome sequencing approach in routine genotyping.

Bayesian Inference Reveals Host-Specific Contributions to the Epidemic Expansion of Influenza A H5N1.

Since its first isolation in 1996 in Guangdong, China, the highly pathogenic avian influenza virus (HPAIV) H5N1 has circulated in avian hosts for almost two decades and spread to more than 60 countries worldwide. The role of different avian hosts and the domestic-wild bird interface has been critical in shaping the complex HPAIV H5N1 disease ecology, but remains difficult to ascertain. To shed light on the large-scale H5N1 transmission patterns and disentangle the contributions of different avian hosts on the tempo and mode of HPAIV H5N1 dispersal, we apply Bayesian evolutionary inference techniques to comprehensive sets of hemagglutinin and neuraminidase gene sequences sampled between 1996 and 2011 throughout Asia and Russia. Our analyses demonstrate that the large-scale H5N1 transmission dynamics are structured according to different avian flyways, and that the incursion of the Central Asian flyway specifically was driven by Anatidae hosts coinciding with rapid rate of spread and an epidemic wavefront acceleration. This also resulted in long-distance dispersal that is likely to be explained by wild bird migration. We identify a significant degree of asymmetry in the large-scale transmission dynamics between Anatidae and Phasianidae, with the latter largely representing poultry as an evolutionary sink. A joint analysis of host dynamics and continuous spatial diffusion demonstrates that the rate of viral dispersal and host diffusivity is significantly higher for Anatidae compared with Phasianidae. These findings complement risk modeling studies and satellite tracking of wild birds in demonstrating a continental-scale structuring into areas of H5N1 persistence that are connected through migratory waterfowl.

Genome-Wide Evolutionary Analyses of G1P[8] Strains Isolated Before and After Rotavirus Vaccine Introduction.

Rotaviruses are the most important etiological agent of acute gastroenteritis in young children worldwide. Among the first countries to introduce rotavirus vaccines into their national immunization programs were Belgium (November 2006) and Australia (July 2007). Surveillance programs in Belgium (since 1999) and Australia (since 1989) offer the opportunity to perform a detailed comparison of rotavirus strains circulating pre- and postvaccine introduction. G1P[8] rotaviruses are the most prominent genotype in humans, and a total of 157 G1P[8] rotaviruses isolated between 1999 and 2011 were selected from Belgium and Australia and their complete genomes were sequenced. Phylogenetic analysis showed evidence of frequent reassortment among Belgian and Australian G1P[8] rotaviruses. Although many different phylogenetic subclusters were present before and after vaccine introduction, some unique clusters were only identified after vaccine introduction, which could be due to natural fluctuation or the first signs of vaccine-driven evolution. The times to the most recent common ancestors for the Belgian and Australian G1P[8] rotaviruses ranged from 1846 to 1955 depending on the gene segment, with VP7 and NSP4 resulting in the most recent estimates. We found no evidence that rotavirus population size was affected after vaccine introduction and only six amino acid sites in VP2, VP3, VP7, and NSP1 were identified to be under positive selective pressure. Continued surveillance of G1P[8] strains is needed to determine long-term effects of vaccine introductions, particularly now rotavirus vaccines are implemented in the national immunization programs of an increasing number of countries worldwide.

Host ecology determines the dispersal patterns of a plant virus.

Since its isolation in 1966 in Kenya, rice yellow mottle virus (RYMV) has been reported throughout Africa resulting in one of the economically most important tropical plant emerging diseases. A thorough understanding of RYMV evolution and dispersal is critical to manage viral spread in tropical areas that heavily rely on agriculture for subsistence. Phylogenetic analyses have suggested a relatively recent expansion, perhaps driven by the intensification of agricultural practices, but this has not yet been examined in a coherent statistical framework. To gain insight into the historical spread of RYMV within Africa rice cultivations, we analyse a dataset of 300 coat protein gene sequences, sampled from East to West Africa over a 46-year period, using Bayesian evolutionary inference. Spatiotemporal reconstructions date the origin of RMYV back to 1852 (1791-1903) and confirm Tanzania as the most likely geographic origin. Following a single long-distance transmission event from East to West Africa, separate viral populations have been maintained for about a century. To identify the factors that shaped the RYMV distribution, we apply a generalised linear model (GLM) extension of discrete phylogenetic diffusion and provide strong support for distances measured on a rice connectivity landscape as the major determinant of RYMV spread. Phylogeographic estimates in continuous space further complement this by demonstrating more pronounced expansion dynamics in West Africa that are consistent with agricultural intensification and extensification. Taken together, our principled phylogeographic inference approach shows for the first time that host ecology dynamics have shaped the historical spread of a plant virus.

Trends and predictors of transmitted drug resistance (TDR) and clusters with TDR in a local Belgian HIV-1 epidemic.

We aimed to study epidemic trends and predictors for transmitted drug resistance (TDR) in our region, its clinical impact and its association with transmission clusters. We included 778 patients from the AIDS Reference Center in Leuven (Belgium) diagnosed from 1998 to 2012. Resistance testing was performed using population-based sequencing and TDR was estimated using the WHO-2009 surveillance list. Phylogenetic analysis was performed using maximum likelihood and Bayesian techniques. The cohort was predominantly Belgian (58.4%), men who have sex with men (MSM) (42.8%), and chronically infected (86.5%). The overall TDR prevalence was 9.6% (95% confidence interval (CI): 7.7-11.9), 6.5% (CI: 5.0-8.5) for nucleoside reverse transcriptase inhibitors (NRTI), 2.2% (CI: 1.4-3.5) for non-NRTI (NNRTI), and 2.2% (CI: 1.4-3.5) for protease inhibitors. A significant parabolic trend of NNRTI-TDR was found (p = 0.019). Factors significantly associated with TDR in univariate analysis were male gender, Belgian origin, MSM, recent infection, transmission clusters and subtype B, while multivariate and Bayesian network analysis singled out subtype B as the most predictive factor of TDR. Subtype B was related with transmission clusters with TDR that included 42.6% of the TDR patients. Thanks to resistance testing, 83% of the patients with TDR who started therapy had undetectable viral load whereas half of the patients would likely have received a suboptimal therapy without this test. In conclusion, TDR remained stable and a NNRTI up-and-down trend was observed. While the presence of clusters with TDR is worrying, we could not identify an independent, non-sequence based predictor for TDR or transmission clusters with TDR that could help with guidelines or public health measures.